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C/^ Anomalous X-ray Pulsars (AXPs) are a class of rare X-ray pulsars whose energy source 

^ has been perplexing for some 20 years^'^'^. Unlike other, better understood X-ray pul- 

y—i sars, AXPs cannot be powered by rotation or by accretion from a binary companion, 

hence the designation "anomalous." AXP rotational and radiative properties are strik- 

^ ingly similar to those of another class of exotic objects, the Soft Gamma Repeaters 

(SGRs). However, the defining property of SGRs, namely their low-energy gamma- 

^ ray and X-ray bursts, have heretofore not been seen in AXPs. SGRs are thought to 

Onbe "magnetars," young neutron stars powered by the decay of an ultra-high magnetic 

S field^'^. The suggestion that AXPs are magnetars has been controversial^. Here we 

d report the discovery, from the direction of AXP IE 1048—5937, of two X-ray bursts 
j> 

that have many properties similar to those of SGR bursts. These events imply a close 

X 

5— I Relationship between AXPs and SGRs, with both being magnetars. 

a 

SGRs are believed to be magnetars because the high magnetic field provides the torque for their 
rapid spin-down, as well as the energy to power their bursts and quiescent X-ray emission^. AXPs 
have been suggested to be magnetars, albeit less active, because of their similar spin periods, rates 
of spin down, location in the Galactic plane, and similar though somewhat softer X-ray spectra to 
those of SGRs in quiesence^. The physical diff'erence between the two classes is unknown, but, in 
the magnetar model, is likely related to the magnitude or distribution of the stellar magnetic field. 
However, the apparent absence of any bursting behavior in AXPs has led to suggestions that they 
could be powered, not by magnetism, but by accretion from a disk of material remaining after the 
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birth supernova event^. If so, the observational similarities between AXPs and SGRs must be purely 
coincidental. 

A program to regularly monitor the AXPs using the Proportional Counter Array (PCA)^ aboard 
NASA's Rossi X-ray Timing Explorer (RXTE) was begun in 1996 in order to determine their long- 
term timing, pulsed flux, and pulse profile stabilities^ As part of this program, motivated 
by the existence of SGR bursts, we also searched the AXP data for bursts (see Fig. 1 caption for 
details) . 

We discovered two highly significant bursts from the direction of AXP IE 1048—5937 in this way. 
The first (hereafter Burst 1) occurred during a 3-ks PCA observation obtained on 2001 October 29 
with chance probability P ~ 6 x 10~^^ after accounting for the number of trials. A second burst 
(hereafter Burst 2) was found in a 3-ks observation obtained on 2001 November 14, with analogous 
probability P ~ 2 x 10~^. No other significant bursts were found toward IE 1048—5937. The total 
PCA time searched for bursts toward this source was 380 ks in observations obtained from 1996-2002. 

The burst profiles are shown in Figure 1. Both arc characterized by fast rises and slow decays 
(see Table 1). Burst 1 appears to have a long, low-level tail that is just above the PCA background 
as determined by intervals selected before and after the bursts (see Fig. 1), while Burst 2 is much 
shorter. Both bursts arrived at the peak of the AXP pulse within uncertainties in burst arrival 
time and definition of pulse peak. The probability of this occurring by random chance is ~1%. We 
note a marginal 2>a) increase in the pulsed flux from IE 1048.1—5937 that commenced with the 
observation in which Burst 1 was detected, and which lasted ~4 weeks. 

To determine the bursts' spectral properties, we first established that neither burst exhibited 
significant spectral evolution with time by computing hardness ratios (the ratio of 10-60 keV counts 
to 2-10 keV counts) for the first 0.5-s and subsequent 1.5-s burst intervals. No significant change in 
hardness was detected, though marginal spectral softening with time was detected after the first 2.5 s 
of Burst 1. Hardness ratios for Burst 1 and 2 for the 1 s after burst onset were 2.8 ±0.8 and 1.0 ±0.3, 
respectively. 

We then fit the spectra from the first 1 s of each burst to two one-component models, a power 
law and a black body (see Table 1). Continuum models provided an adequate characterization of the 
Burst 2 spectrum but not of the Burst 1 spectrum. As seen in Figure 2, the spectrum for the 1 s 
after the Burst 1 onset exhibits a feature near 14 keV. This feature is clear in all binning schemes and 
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is prominent throughout the first ~1 s of the burst. No known PCA instrumental effect produces a 
feature at this energy (K. Jahoda, personal communication). 

Due to the wide (~1°) field-of-view (FOV) and lack of imaging capabilities of the PCA, we 
cannot verify that the bursts originated from the location of the AXP. The low peak X-ray fluxes 
of the events (see Table 1) preclude determining the source's location using data from other, better 
imaging instruments that were contemporaneously observing the X-ray sky, such as the RXTE All 
Sky Monitor, or the Wide Field Camera aboard BeppoSAX. We must therefore consider other possible 
origins from the bursts before concluding they were from the AXP. 

The bursts' short rise times (Table 1) require emission regions of less than a few thousand km, 
implying a compact object origin. So-called Type I X-ray bursts are a well-studied phenomenon that 
result from unstable helium burning just below the surface of a weakly magnetized neutron star that 
is accreting material in a low-mass X-ray binary (LMXB)^^. However, Type I bursts from an LMXB 
in the same FOV as IE 1048.1—5937 are unlikely to explain our observed bursts because (i) the burst 
rise times arc much shorter than those of Type I bursts; (ii) the burst spectra are much harder than 
those of Type I bursts; (iii) Burst 2 shows no evidence for spectral softening with time and no Type 
I burst has ever exhibited a spectral feature like the one detected in Burst 1; (iv) the bursts are 
extremely faint, implying a source location well outside the Milky Way for Type I burst luminosities 
(v) there are no known LMXBs in the FOV^^. Type II X-ray bursts^^ are a much rarer and less well 
understood phenomenon observed thus far in only two sources, both accreting binaries. The bursts 
we have observed are unlikely to be Type II bursts from an unknown X-ray binary in the PCA FOV 
because (i) of the rarity of such events; (ii) Type II bursts have longer rise times than do our bursts; 
(iii) no Type II burst has exhibited a spectral feature like that seen in Burst 1. 

Classical gamma-ray bursts (GRBs) sometimes exhibit prompt X-ray emission that can have 
temporal and spectral signatures similar to those wc have obscrvcd^^. However, the likelihood of two 
GRBs occurring within 1° of each other is small, and GRBs are not known to repeat. Conservatively 
assuming GRB spectral model parameters that result in low gamma-ray fluxes and extrapolating 
the GRB rate^^ as measured with the Burst and Transient Source Experiment (BATSE)^^ assuming 
homogeneity below the BATSE threshold, we estimate a probability that these events are unrelated 
GRBs that occured by chance in the same RXTE FOV during our IE 1048.1—5937 monitoring 
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observations (conservatively neglecting that they occured within two weeks of each other) of ^ 9 x 
10-^ 

The observed burst properties are in many ways similar to those seen from SGRs"*^^ . The fast rise 
and slow decay profiles are consistent with SGR time histories, as are the burst durations (neglecting 
the long, low-level tail of Burst 1). Both AXP and SGR bursts are spectrally much harder than is 
their quiescent pulsed emission. The burst peak fluxes and fluences fall within the range seen for 
SGRs, and the spectrum of Burst 2 is consistent with SGR burst spectra of comparable fluence. 
Burst 1 has characteristics unlike nearly all SGR bursts, specifically its long tail and spectral feature. 
However, we note that a single event from SGR 1900+14 was shown^^'^*^ to possess each of these 
properties. The marginal increase in the pulsed fraction that we observed at the burst epochs is 
consistent with SGR pulsed flux increases seen during bursting episodes^^. Finally, the fact that in 
spite of several years of monitoring, the only two bursts detected occurred within two weeks of each 
other suggests episodic bursting activity, the hallmark of SGRs. Thus, the characteristics of these 
events match the burst properties of SGRs far better than any other known burst phenomenon. 

In the magnetar model for SGRs^'^, bursts are a result of sudden crustal yields due to stress from 
the outward diffusion of the huge internal magnetic field. Such yields cause crust shears which twist 
the external magnetic field, releasing energy. Thompson &; Duncan^ who, upon suggesting that AXPs 
are also magnetars, predicted X-ray bursts should eventually be seen from them. By contrast, in no 
AXP accretion scenario, whether binary or isolated fall-back disk, are SGR-like bursts expected. 

The large 14-keV spectral line in Burst 1 is intriguing. An electron cyclotron feature at this 
energy E implies a magnetic field oi B = 2TrmcE/he 2± 1.2 x 10-*^^ G (where m is the electron mass, 
c is the speed of light, h is Planck's constant, and e is the electron charge), while a proton cyclotron 
feature implies B ~ 2.4 x 10^^ G. The former is significantly lower than is implied from the source's 
spin-down and is typical of conventional young neutron stars, rather than magnetars. The latter is 
higher than is implied by the spin-down yet reasonable for the magnetar model as the spin-down 
torque is sensitive only to the dipolar component of the magnetic field. 

Why do the burst rates of AXPs and SGRs differ so markedly, in spite of their common magnetar 
nature? One possibility is that AXP internal magnetic fields are much larger than those of SGRs; 
if so, AXP crusts can undergo plastic deformation rather than brittle fracturing''. However, this is 
opposite to what is inferred from the two classes' spin-down rates, suggesting the latter is an unreliable 
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internal field indicator. This could help reconcile the contrasting radiative properties of AXPs and 
apparently high-magnetic field radio pulsars^^. It also suggests that AXPs are SGR progenitors, with 
bursting behavior commencing as the field decays. This is consistent with the smaller AXP ages 
implied by their more numerous associations with supernova remnants^^, but does not explain why 
AXPs and SGRs have similar spin period distributions, since AXPs spin down as they age^^. This 
aspect of magnetar physics remains a puzzle. 

Received 2 February 2008; Accepted draft. 
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Table 1 



AXP Burst Timing and Spectral Properties 



Burst 1 



Burst 2 



Burst day, (MJD) 
Burst start time, (fraction of day, UT) 
Burst rise time, tr (ms) 
Burst duration, Tgo (s) 
Burst phase 



Temporal Properties 

52211 



0.2301949(24) 



2ltl 



51 



+28 
-19 



52227 
0.836323379(68) 
5.91?:° 

9 n+4-9 
^•'-'-0.7 



-0.018 ± 0.034 0.051 ± 0.032 



Fluxes and Fluences 

485 ±118 



Tgo fluence (counts) 
Tgo fluence (xl0~^° erg cm~^) 
1-s fluence (counts) 
1-s fluence (xlO~^'^ erg cm~^) 
Peak flux for 64 ms (xlO~^° erg s~^ cm ^) 
Peak flux for tr ms (xlO"^*^ erg cm ^) 

Spectral Properties 

Power law: 

power law index 

power law flux (xlO~^° erg s"^ cm~^) 
line energy (keV) 
line width, a (keV) 
line flux (xlO"^'' erg s~^ cm~^) 
reduced /degrees of freedom 
Black body: 
kT (keV) 

black body flux (xlO"^'' erg s^^ cm~^) 
line energy (keV) 
line width (keV) 
line flux (xlO"^'^ erg s"^ cm~^) 
reduced /degrees of freedom 



20.3 ± 4.8 
117 ± 13 



0-89lJ:?i 
13.9 ± 0.9 



2.2 



-1.3 
-1.0 



1.24/15 



3.9ii:? 

2.411? 



14.2 



+1.1 

-1.2 



2 1+1-5 



3-71?:^ 

1.23/15 



101 ± 15 
5.3 ± 1.2 

69 ± 10 

A n+3.5 
^■^-0.8 

26lf 



1 qo-l-0.75 

l.OO_0,g2 

4.olg:| 



0.77/5 



3.6lf;2 
3.810 



1.66/5 
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Table Caption 

Table 1. The uncertainty on the burst start time is the burst rise time tr and is given in 
parenthesis as the uncertainty in the last digits shown. The burst rise times were determined by 
a maximum hkehhood fit to the unbinned data using a piecewise function having a hnear rise and 
exponential decay. The burst duration, T90, is the interval between when 5% and 95% of the total 
2-20 keV burst fluence was received. The background regions used for both the duration and spectral 
analyses are shown in Figure 1. Burst phase is defined such that the peak of the periodic pulsation 
is at phase 0/1. All fiuences and fluxes are in the 2-20 keV range. Tgo fluences in cgs units were 
calculated assuming a power-law spectral model and spectral grouping that demanded a minimum 
of 20 counts per spectral bin. The 1-s fluences in cgs units correspond to the fluxes found in the 
spectral modeling. The spectral rcbinning method used in all spectral modeling for Burst 1 was to 
group the 256 PCA channels by a factor of 4, while for Burst 2, we demanded at least 20 counts 
per spectral bin. Peak fluxes on the short time scales were determined by scaling the 1-s fluxes by 
number of counts. For all spectral fits, the equivalent neutral hydrogen column density was held 
fixed at 1.2 x 10^^ cm~^, the value determined from recent XMM observations^^. Spectral fits were 
determined for the first 1 s after the burst start times. Spectral modeling was done using photons in 
the 2-40 keV range. Response matrices were created using the FTOOL pcarsp ^. Uncertainties in the 
Table are 68% confidence intervals, except for those reported for the cgs-unit fluences and fluxes, as 
well as the spectral model parameters, for which we report 90% confidence intervals. 



^http://heasarc. gsfc.iiasa.gov/docs/xto/recipes/pca_response. html 
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Burst 1 Burst 2 




Time (seconds from onset) 



Time (seconds from onset) 
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Figure 1 Caption 

Lightcurves for the observed bursts. The RXTE AXP data set consists of short (~3 ks) snapshots, 
as well as longer archival observations, all taken in the PCA GoodXenonwithPropane mode, which 
records photon arrival times with l-/is resolution, and bins photon energies into 256 channels. Time 
series were initially created with 31.25-ms resolution from photons having energies in the range 2— 
20 keV for each PCA Proportional Counter Unit (PCU) separately, using all xenon layers. Photon 
arrival times at each epoch were adjusted to the solar system barycentre. The resulting time series 
were searched for significant excursions from the mean count rate by comparing each time bin value 
with a windowed 7-s running mean. Bursts were identified assuming Poissonian statistics, and by 
combining probabilities from the separate PCUs. Left Panel: Background subtracted 2-20 keV light 
curves for Burst 1, displayed with 2-s time resolution. The solid horizontal lines before and after the 
bursts are the boundaries of the pre- and post-background intervals used for calculating Tgo and for 
spectral modeling. The Tgo interval is shown as a horizontal dashed line. Right Panel: Same but for 
Burst 2, and with 0.5-s time resolution in the main panel. The insets show the peak of each burst 
with 31.25-ms time resolution. We verified that there was no significantly enhanced signal from PCA 
events not flagged as "good" at the times of the bursts (such as those that do not enter through the 
PCA aperture) in the RXTE "Standard 1" event files. We also verified that both events were clearly 
detected in all operational PCUs. Hence the events are unlikely to be instrumental in origin. 
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Figure 2 Caption 

X-ray spectrum in the 2-40 keV range for the 1 s after the onset of Burst 1. For all spectral analyses, 
we first created spectral files having 256 channels across the full PCA energy range (~0.2-60 keV), 
although subsequent factor of 4 rebinning was necessary because of the paucity of burst counts. The 
burst and background intervals were used as input to the X-ray spectral fitting package XSPEC^^ 
vl 1.1.0^. The spectrum of the first 1 s after Burst 1 onset is not well characterized by any continuum 
model. The best fit power-law plus Gaussian line model is shown as a solid line. The F-test shows 
that the addition of a line of arbitrary energy, width and normalization to a simple power law model 
improves the fit significantly, with a chance probability of this occuring of 0.0032. Monte Carlo 
simulations in XSPEC were done to verify this conclusion: 10,000 simulations of similar data sets were 
produced assuming a simple power-law energy distribution, then fit with a power law plus Gaussian 
line of arbitrary energy, width and normalization. This procedure is conservative, since it ignores that 
the observed large line has flux comparable to the measured continuum. In 10,000 trials, we found 1 
trial with the same or smaller chance occurence probability as judged by the F-test, indicating that 
the probability of the line we observed being due to random chance is < 0.0001. We repeated this 
procedure for data having a black-body spectrum, with similar results, namely the probability of the 
line being due to chance is < 0.0008. The spectrum also shows possible additional features at ~7 keV 
and ~30 keV (suggestive of of lines at multiples of 1, 2 and 4 of ~7 keV). These additional features 
are not apparent in all binning schemes and are not statistically significant. 



^http: / /xspec. gsfc.nasa.gov 



